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ABSTRACT 

We study moderate gravitational lensing where a background galaxy is magnified substan- 
tially, but not multiply imaged, by an intervening galaxy. We focus on the case where both the 
lens and source are elliptical galaxies. The signatures of moderate lensing include isophotal 
distortions and systematic shifts in the fundamental plane and Kormendy relation, which can 
potentially be used to statistically determine the galaxy mass profiles. These effects are illus- 
trated using Monte Carlo simulations of galaxy pairs where the foreground galaxy is modelled 
as a singular isothermal sphere model and observational parameters appropriate for the Large 
Synoptic Survey Telescope (LSST). The range in radius probed by moderate lensing will be 
larger than that by strong lensing, and is in the interesting regime where the density slope may 
be changing. 

Key words: Gravitational lensing - galaxies: elliptical and lenticular - dark matter - Galaxies: 
structure, Galaxies: formation 



1 INTRODUCTION 

Gravitational lensing is usually divided into microlensing, strong 
lensing, and weak gravitational lensing. Microlensing refers to the 
temporal change of magnification of a background source lensed by 
an intervening object (usually stars). Strong lensing occurs when 
multiple images or strong distortions (e.g. giant arcs) of a back- 
ground source are caused by an intervening galaxy or cluster. For 
weak lensing the distortion of a background source by the lensing 
object is much more subtle. For example, a circular source will be 
lensed into an ellipse with ellipticity of a few percent. Such sub- 
tle deviations have to be inferred statistically by averaging over a 
large number of background galaxies. All these fields in gravita- 
tional lensing have been extensively studied, with diverse applica- 
tions ranging from cosmology to the detection of extrasolar planets 
(for a review on these topics, see Schneider et al. 2006 and refer- 
ences therein). 

In this work, we shall explore the intermediate regime which 
we term as "moderate gravitational lensing". In this case, the mag- 
nification is still significant, but no multiple images occur. In the 
context of clusters of galaxies, Williams & Lewis ( 1998jl and |Futa-] 
|mase et aL] fl998) have considered a class of gravitationally lensed, 
highly magnified, yet morphologically regular images (originally 
motivated by observations of the object cB58,|Yee et al.|1996}|Seitz| 
|et al.|1 998). More recently, Sonnenfeld et al. (20TTi discussed the 
magnification effect in with a focus on how to remove the mass- 
sheet degeneracy. In this paper, we shall explore the case where 



an elliptical galaxy is lensed by a foreground elliptical galaxy. 
Many cases are expected to be found in future surveys (e.g. by 
PANSTARRjQ and LSSlQl where hundreds of millions of galax- 
ies will be imaged. Many pairs of galaxies that are close to each 
other but at different redshifts will be discovered. Most of these will 
not be multiply-imaged. At very large separations, a weak-lensing 
galaxy-galaxy analysis will be appropriate; at very small separa- 
tions, multiple images may form. In this paper we carry out Monte 
Carlo simulations of moderate galaxy-galaxy lensing at intermedi- 
ate separations to illustrate the signatures, and what we can learn 
from these. 

The paper is organised as follows. In Section [2] we describe 
the simple singular isothermal sphere model we use. In Section 
[3] we present the main results, including the predictions for the 
optical depth, isophotal distortions, and systematic offsets in the 
fundamental plane and Kormendy relation. We present a summary 
and discussion in Section]?] Throughout this paper, we adopt a flat 
ACDM cosmology model with fl m> o = 0.27 and fl A ,o = 0.73 for 
the matter and cosmological constant, and the Hubble constant is 
written as H = lOOh km s -1 Mpc~' with h = 0.705 (Komatsu et 
|al.|2009) . 
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2 THE LENS AND SOURCE MODEL 

In this section, we describe the lens mass model, the source pop- 
ulation and the surface brightness models of foreground and back- 
ground galaxies which we will use for Monte Carlo simulations of 
galaxy pair catalogues. 



galaxies at redshift 2. Below this truncation redshift, we assume a 
constant comoving number density of lenses, implying the physi- 
cal volume density evolves as OoxCzrf, cr) = ©gj,(c)(1 + z) 3 , where 
the factor (1 + z) 3 accounts for the expansion of the universe. Since 
most of the foreground (lensing) galaxies are below redshift 1, this 
assumption will not have much effects on our results. 



2.1 Lens mass model 

For simplicity, we model the mass profile of the lensing galaxy as 
a singular isothermal sphere (SIS). This simple model is analyti- 
cally tractable, and fits the observational data well (e.g., Koopmans 
|2009] >. The surface mass density can be modelled as 

K?) = ~. (1) 

where f is the (physical) distance in the lens plane, and cr is the 
one-dimensional velocity dispersion. The angular Einstein radius 
of an SIS is given by 
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where c is the speed of light, £> ds is the angular diameter distance 
between the lens and the source, and D s is the angular diameter 
distance to the source. 

A point source at unlensed angular position < 9 S < 9 E 
from the galaxy produces two images at the angular positions 
#i,2 = 9 S ± 9e- The magnification of each image and the total abso- 
lute magnification are given by 
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Notice that the total magnification for the two images exceeds two. 

A point source at angular separation # s > 9 E produces only 
one image at the angular position 6 = 9 S + 9 E , with magnification 
(#e+#s)/&,, giving a value between 1 and 2. Thus if a source is a few 
Einstein radii away from the line of sight, it may still experience 
substantial magnification and differential magnification, and have 
observable signatures. 



2.2 Lensing probability 

For a background elliptical galaxy at redshift z s , the lensing proba- 
bility (optical depth) can be obtained as follows 
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where <x cl . is the cross-section for moderate lensing, t is the cosmic 
time, and QAzi, cr) is the lens velocity dispersion function (e.g. 
Turner, Ostriker & Gott 1984, Gott, Park & Lee 1989, Fukugita, 
Futamase & Kasai 1990). 

The lens velocity dispersion function is modelled by a modi- 
fied Schechter function (Ch oi et al.|2007) : 
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where®* = 8.0xl(r 3 h 3 Mpc~ 3 , a a = 2.32±0.10, y S IT = 2.67±0.07, 
and cr* = 161 ± 5 kms~' . We limit the velocity dispersion to the in- 
terval 70 ~ 400 kms~' l |Loeb & Pe ebles 2003). Most massive early- 
type galaxies were already assembled at z < 1 ; beyond this redshift, 
the density of early-type galaxies declines significantly (e.g., |RerT| 
|zini|2006) , and so we truncate the redshift of early-type background 



2.3 Elliptical galaxy luminosity function 

The luminosity function (LF) of early-type galaxies is modelled by 
the |Schechtef|(T976l > form 

L (L).L = O^A) l \xp(-A)g. (6) 



We adopt values from |Choi et aT] ([2007 1 that were obtained us- 
ing a galaxy sub-sample of the SDSS Data Release 5 with M* - 
5 log 10 h = -20.23 ± 0.04, a L = -0.527 ± 0.043, and <££ = 
0.71 x 10~ 2 fr~ 3 Mpc~ 3 , appropriate for the SDSS r-band. 

It is more convenient for our calculation here to use the lumi- 
nosity function with absolute magnitude M rather than L. Eq. |6]( 
can be rewritten as 



0> L (M)dM = 0.4 In 10x0^10 
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where M 4 is given below eq. |6](. 

Each survey has its own flux limit which corresponds to a 
magnitude limit in a photometric band (e.g. the SDSS r band), de- 
noted by m r . The absolute limiting magnitude, M r , of a galaxy at 
redshift z and galactic coordinates (/, b), can be constructed from 
the apparent magnitude limit m r as follows: 



M r = m r - DM(z) - K r (z), 



(8) 



where DM(z) is the distance modulus, and K,-{z) is the K-correction 
( |Hogg|1999'l >. Here we use the K r {z) given by |Choi et aL] ( |2007) , 
and we have ignored dust extinction by the Milky Way and the 
foreground galaxy. 

The lensing probability averaged over all source redshifts is 
given by (e.g., |Wyithe et al.|2010) 

Nq 
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where A'g is the total number of background elliptical galaxies, Vc 
is the comoving volume, and 'i'o.s (M r (z s ), z s ) is the comoving num- 
ber density of elliptical galaxies brighter than M r (z s ) at redshift z s . 
N G can be calculated as 
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andT G ^(M r (z s ),z s ) as 

-MM) 



Va,s(M r (z 



Oi 



(M)dM, 



(10) 



(11) 



where <b E (M) is the background elliptical galaxy luminosity func- 
tion (see eq.|7J, and M r (z s ) is the absolute limiting magnitude M r at 
redshift z s corresponding to the magnitude limit m r . To be conser- 
vative and for simplicity, we use the LSST single-visit depth limit 
m r = 24.7 in our optical depth calculation. 

2.4 Lens and source surface brightness profiles 

We model the background source (and the lens) as an elliptical 
galaxy with a de Vaucouleurs surface brightness profile: 
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where R is the two-dimensional radius, R c is the effective radius 
within which half of the light is enclosed and / e is the surface 
brightness at the effective radius. We quote the effective radius as 
the geometric mean of the major and minor axes, R e = -\/a e b e . The 
total apparent magnitude mj, average magnitude within effective 
radius (p) t and effective radius R c can be related by i Scod eggio et| 
|al.|1998| > 



m T = -2.51og,o(L T ) - 48.6 = </z> e - 1.9954 - 51ogi (Re), 



(13) 



where L T is the total flux of the elliptical galaxy and R e is measured 



2.5 Simulating images of early-type galaxies 

In an image of an early-type galaxy, in addition to the signal of 
the source, we also have the Poisson noise in the sky background, 
readout noise of the CCDs and the dark current. The signal-to-noise 
ratio is given by (McM aster et al.|2008"l > 

S_ = N s t 

N 



tJn s i + n(N Ay t + N x t + Nl) 



(14) 



where N s is the number of photoelectrons from the source per unit 
time, N^ y is the intensity of sky background in photoelectrons per 
pixel per unit time, Noc is the dark current in electrons per pixel 
per unit time, N R is the readout noise in electrons per pixel, n is the 
number of pixels covered by an image of an early-type galaxy and 
t is the exposure time. 

To proceed further, we need to adopt a concrete case for the 
parameters in eq. ( |14) , Here we restrict ourselves to LSST The tele- 
scope will be sited in Cerro Pachon, Chile, with excellent median 
free-air seeing of 0."7 in the r band. Correspondingly, the pixel 
scale is chosen to be 0."2. The telescope will repeatedly survey the 
sky covering an area of 20000 square degrees with a cadence of 
~ 7 days. The integration time (t) for each individual frame is 15 
seconds. We assume that the readout noise is the same for all expo- 
sures. For the Poisson noise in the sky background, readout noise 
of the CCDs and the dark current, we use the mean values given 
by the LSST Exposure Time Calculate^] Using these numbers, we 
verified that the photon numbers we obtain differ from those given 
by the LSST Exposure Time Calculator by less than 3%. 

Fig.[T] shows an example of a simulated moderately lensed 
background elliptical galaxy over-plotted with the best-fit ellipses 
from the IRAF task ELLIPSE (for more details, see j |3.2[ (. Notice 
that we have not convolved the image with seeing. 

To simulate a population of early-type galaxies, we use the 
method given in Appendix A of Ber nardi et al.| ( |2003a) . Briefly, 
each galaxy is randomly assigned an absolute magnitude, M r , ef- 
fective radius R e and velocity dispersion a (V in the notation of 
Bernardi et al. 2003a), appropriate at redshift z = 0. This proce- 
dure ensures the generated galaxies satisfy the observed correla- 
tions among observed properties, including the fundamental plane 
as given by |Bernardi eULl < |2003b| 



log 10 o- + 0.2«ju>e- 20.09). 



(15) 



We assign a typical error of 0.03 dex in log <r and 0.02 dex in 
log 10 i? c . The scatter in the fundamental plane is quite small, on 
the order of as 0.08, and may be mostly due to measurement errors. 



http://dls. physics. ucdavis.edu:8080/etc4_3work/servlets/LsstEtc. html 
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Figure 1. A simulated lensed image of a background elliptical galaxy, with 
a scale of 0."2 per pixel. The foreground lens is at redshift 0.3 and the 
background elliptical galaxy (on the right) is at redshift 0.6. The size of the 
image is 30. "0 by 25. "6. The distance between the two galaxy centres is 
8.2 arcsec. The ellipses are obtained by running the IRAF task ELLIPSE 
(see j |3.2[ . Notice the small changes in the position angle and twists in the 
isophotes. 



The Kormendy ( 1977) relation between the surface brightness 
and effective radius is given by ( jBernardi et al. 2003a as quoted in 
IQohama et al.|2009| > 



<^> e = 2.04 log 10 tf c + 18.7, 



(16) 



where (yu) c is the mean surface brightness (in mag/arcsec 2 ) within 
the effective radius (R e , in units of kpc). The scatter in this relation 
is substantially larger than that in the fundamental plane (see §3.4) . 

For LSST, the mean redshift of the galaxy sample is about 0.8 
(see the LSST science book for details). However, to be specific and 
somewhat conservative, we use a lens redshift of 0.3 and a source 
redshift of 0.6. In reality, galaxy pairs will have a redshift distribu- 
tion, but we can in principle bin the data and study the quantities 
we are interested, and so this simplification does not change the re- 
sults of the paper. Notice that for the optical depth calculation, we 
still assume the lenses and sources have a redshift distribution, as 
described in ^2] 

For each random galaxy generated for redshift 0, we shift it 
to redshift z by taking into account its evolution |Bernardi et al.| 
|2003a) 



M r (z) = M r (z = 0) - Qz, 



(17) 



where Q = 0.85 for the r band. 

We assume a total exposure time of 6000 seconds (i.e., 400 
exposures of 15 seconds). We further select only galaxies which 
satisfy R e > 0."7, as a way to crudely account for the effects of 
seeing. In total, we generate 11623 pairs of foreground and back- 
ground galaxies which satisfy these conditions. The ellipticity of 
the background galaxy is randomly drawn from the axis ratio distri- 
bution as given in |Choi et al.| ( |2007] > (see their Fig. 13). The lensed 
background galaxy images are obtained using the lens equation. 
The following analyses are based on this sample, although we ap- 
ply a further cut in the isophote shape analysis (see below). 
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Figure 2. The optical depth, t, as a function of the redshift, z s . The cross- 
section for moderate lensing is taken as the area between 16g to 10$e- 



Figure 3. The probability density distribution of the effective radii of the 
background elliptical galaxies, all at redshift 0.6. At this redshift, 1 arcsec 
corresponds to 4.7 h~ i kpc. 



3 RESULTS 
3.1 Optical depth 

In Fig. [2] we present the optical depth as a function of redshift. In 
this exercise, we take the moderate lensing cross-section for each 
galaxy as 9e to 10#e (corresponding tOyU = 2 and 1.1). Not surpris- 
ingly the optical depth increases with increasing redshift, reflecting 
the fact that the number of intervening galaxies increases for more 
distant background sources. The optical depth is about 0.0029 at 
z = 0.5, 0.016 at z = 1.0, and 0.062 at z = 2.0. The optical depth 
of averaged over all the background sources is about 0.0254. This 
probability will scale linearly with the cross-section we adopt. For 
example, if the we take the cross-section from 29 E to 56 E (corre- 
sponding to fj = 1.5 and 1.2) then the probability will be reduced 
by a factor of roughly 5 to about 0.005. The precise choice of cross- 
section will depend on the depth, seeing etc. of observations. 

We have also calculated the magnification bias (Turn eFet al.| 
|1984| >, which can be large for multiply-imaged quasars or galaxies. 
In our case, we find the magnification bias to be modest, only about 
1.16. The magnification bias is small because the magnification in 
our case is by definition modest. 

As we have mentioned in j |2.5| the redshift of the background 
sources is truncated at z = 2.0. With this limitation, the surface 
number density of early-type galaxies is about 10 per square ar- 
cmin. For LSST, it will carry out a survey of 20000 deg 2 of the sky, 
and about 700 million elliptical galaxies will be detected. In Fig. [3] 
we present the distribution of the effective radius of the background 
elliptical galaxy sources. The number of elliptical galaxies drops 
steeply for larger effective radii. About 54.1% of the background 
sources have effective radii larger than 0."7, the median free-air 
seeing at the site of LSST. Under this criterion, we estimate that 
there are about 9.5 million moderate galaxy lensing cases that will 
be observed by LSST to the single- visit depth (m,- = 24.7 mag). 

The predicted angular Einstein radius is shown in Fig. [4] As 
can be seen, the distribution peaks around 0.8 arcsec and has an 
extended tail out to larger separations. The distribution is similar to 
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Figure 4. The probability density distribution of the angular Einstein radius 
&E of the simulated lensing cases. 



the observed one from the CLASS survey, even though the source 
and lens distributions there are somewhat different (see Fig. 1 1 in 
|Browne et al.|2003) . 



3.2 Isophotal distortions 

Surface photometry is performed on the simulated lensed images 
of the background elliptical galaxy by utilizing the IRAF task 
ELLIPSE l |Jedrzejewski|1987) . When running the ELLIPSE on the 
lensed images, the geometric centre, elllipticy and position angle 
are all allowed to vary freely, with a logarithmic step of 0.05 in the 
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Figure 5. Coefficients of the Fourier series a^ja, 04/0, b}/a and /34/a ob- 
tained by the IRAF task ELLIPSE for the galaxy shown in Fig.[T] The hori- 
zontal line in each panel indicates the value for a perfect ellipse. 



semi-major axis. The isophotes of the lensed images can be mea- 
sured in different ways. Here we present the output of the task EL- 
LIPSE. A detailed description of these measurements can be found 
in |Hao et al.| l |2006) , we only repeat the essentials here. The inten- 
sity along the ellipse is expanded in Fourier series 



1(0) = To + y^CAi cos n9 + B n sin nff), 



(18) 



where I is the intensity averaged over the ellipse, and A„ and B„ 
are the higher order Fourier coefficients. If an isophote is a perfect 
ellipse, then all the coefficients, (A n , B„), n = 1, • • •, 00 will be ex- 
actly zero. We will use two quantities a„ and b„ extensively, which 
are related to the ELLIPSE outputs A„ and B„ by 

i = da, (19) 

a ya 

where y = dl/dR is the local radial intensity gradient, and a is the 
semi-major axis length. 

For the galaxy shown in Fig.[T] the coefficients of the Fourier 
series a 3 /a, a 4 /a, b^/a and b^/a as a function of the semi-major 
axis radius are shown in Fig.|5] The horizontal line in each panel is 
the value for a perfect ellipse. Notice that both a^/a and b^/a show 
significant deviations from the isophotes of a perfect ellipse. 

For the same galaxy, we present the distribution of the appar- 
ent isophotal axis ratio q, position angle ip in Fig.[6]. The horizontal 
line presents the ellipticity (position angle) before lensing. There 
are some small but systematic changes in these two quantities of 
a few percent for the example shown here. The 04 /a and £14/12 pa- 
rameters describe whether an elliptical galaxy is disky or boxy, but 
we did not incorporate these into our image simulations. Instead we 
choose to focus on the a^/a and b^/a parameters. 

In Fig|7] we present the distribution of the coefficients of the 
Fourier series a^ja, 04/0, b^/a and b 4 /a for cases where the sepa- 
rations between the foreground galaxy and background galaxy are 
between 9 E + R e s » 2."0 and 10. "0, where R e s is the effective ra- 
dius for the source. For each galaxy, these coefficients are the mean 
value within the region between 2r s and l.SRsa, where r s a 0."7 is 
the median seeing radius and R 50 is the Petrosian half-light radius. 
For SDSS, R50 = 0.7 IR C for early-type galaxies described by the 
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Figure 6. The probability density distributions of the apparent isophotal 
axis ratio (q), position angle (ip, in units of degrees) for the galaxy shown in 
Fig-Q] 



R 1 ' 4 law ( |Hao et al.|2006[|Graham et al.|2005| . In the same plot, 
these distributions for separations between 2."0 and 6."0 and be- 
tween 2."0 and 10. "0 are also shown. As expected, for background 
sources with smaller separations from the lens centre, the distribu- 
tions of the coefficients are broader because the lensing effects are 
stronger in such cases. Notice that in the above analysis, we require 
2r s > 1.5^50, this condition reduces the pair of galaxies from 1 1623 
to 9094. 

Both the a 3 /a and b 3 /a distributions show large deviations 
from those of normal elliptical galaxies as studied by |Hao et al.| 
(20061. This is further illustrated in Fig. [8] here we plot the two- 
dimensional distributions of these two quantities rather than as his- 
tograms from projection. We can see that the unlensed galaxies are 
mostly at the centre, while the lensed galaxies clearly show much 
larger scatters than the unlensed ones. Thus they can be used as 
effective indicators whether moderate gravitational lensing has oc- 
curred in a pair. 

3.3 Fundamental plane 

Gravitational lensing preserves the surface brightness of the source 
because of Liouville's theorem, but changes the apparent solid an- 
gle of a source. Following this, a natural expectation is that the fun- 
damental plane of background elliptical sources should be affected 
by moderate gravitational lensing. 

To obtain the lensed and unlensed effective radii and surface 
brightness, R cE and R c , we fit the surface brightness profile of the 
lensed images of background elliptical galaxies with the R [/4 law. 
In this model, we fix the centre of the galaxy, after which there are 
four remaining free parameters, j/ c , R c , tp, q}. The x 1 function is 
given by 



X z (Ie,Rc,<P,q) = 2_j t— 

i=\ °i,0 



(20) 



where N p is the number of pixels of the lensed image. J, o is the 
photon number in the z' th pixel, <x, is the Poisson noise cr, = 



I la- The degree of freedom is iVdof = iV p - 4. We use the rou- 
tine gsl_multimin_fminimizer_nmsimplex in the GNU Scientific Li- 
brarjQto find the^ 2 minimum and obtain the effective radius; the 
mean surface brightness is then obtained by averaging the surface 
brightness within R s . 

Fig.|9] shows the distribution of differences in effective radii 



http://www.gnu.org/software/gsl/ 
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Figure 7. The probability density distributions of the coefficients of the 
Fourier series a^/a, a^fa, b}/a and b^/a obtained by the IRAF task 
ELLIPSE. The black solid histogram in each panel shows the distribution 
of the coefficients for a nearby sample of early-type galaxies from |Hao et| 
| al. 1(2006) . The green line is for moderately lensed galaxy pairs within sep- 
aration between 2 to 6 arcsec while the red solid line is for galaxies from 2 
to 10 arcsec. 
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Figure 8. The scattered plot for the simulated galaxy pairs in the plane 03/0 
vs. i>3 /a. The red dots are for unlensed galaxies while the black ones are for 
lensed ones. In total, there are 9094 pairs of galaxies. 



between the lensed and unlensed sources, denned as Alog 10 yi? e = 
log 10 R c L - log 10 R c . The mean and median values of Alog[ R c are 
0.03 and 0.03 respectively. These values will depend on the separa- 
tion range we take for the galaxy pair, here we have selected pairs 
with separations between ~ 2 to 10 arcsec. 

Similarly, Fig.fTU] shows the distribution of differences in the 
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Figure 9. The probability density distribution of differences in effective 
radii between the lensed and unlensed sources in our simulated sample. 
Notice the systematic shift to a larger effective radius. 
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Figure 10. The probability density distribution of differences in mean sur- 
face brightness within effective radius between the lensed and unlensed 
sources in our simulated sample. 



mean surface brightnesses within the effective radius between the 
lensed and unlensed sources, defined as A(yu) e = </j) c .l - 
where <yu) c _L and </i> e are the mean surface brightness within the 
lensed and unlensed effective radius of the lensed images. The 
mean and median differences are -0.0003 and -0.0003, which are 
much smaller than those for the effective radius, i.e., gravitational 
lensing does not affect the mean surface brightness within the ef- 
fective radius. 

The fundamental plane of the lensed elliptical galaxies is pre- 
sented in Fig.[TT| The black solid line shows the fundamental plane 
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Figure 11. Fundamental plane of the lensed elliptical galaxies. The black 
solid line shows the fundamental plane for the background elliptical galax- 
ies. The red dots are for the lensed elliptical galaxies with the best regression 
line shows as the red dashed line. The lens galaxies are at redshift z.\ = 0.3 
and the source galaxies at z s = 0.6. 

used to produce the background elliptical galaxies (black dots) 
while the red symbols present the lensed elliptical galaxies. The 
best regression line for the lensed sample is shown as the dashed 
red line. On average an upward offset of about 0.03 is introduced 
by gravitational lensing. 

3.4 Kormendy relation 

To derive the systematic shifts in the fundamental plane, in princi- 
ple we need to measure the velocity dispersion from spectroscopic 
observations. This will be very time consuming given the large 
number of early-type galaxies available from LSST. However, it 
is easier to measure the Kormendy ( 1977 1 relation between the sur- 
face brightness and effective radius, which can be obtained from 
imaging data and an approximate photometric redshift. In Fig.|12| 
we show the Kormendy relation, of the lensed and unlensed ellip- 
tical galaxies. On average, an upward offset of about 0.03 in the 
Kormendy Relation is seen in Fig. [TT| As might be expected, this 
shift is similar to that in the fundamental plane. This offset is much 
smaller than the scatter orthogonal to the Kormendy relation (about 
0.166) and the scatter in the direction of \ogR c (about 0.182). 



4 CONCLUSION AND DISCUSSION 

In this paper we have investigated the moderate lensing of back- 
ground elliptical galaxies by intervening elliptical galaxies. We find 
the the optical depth for moderate lensing is on the order of ~ 1%. 
We have also performed image simulations based on the design 
specifications of LSST and obtained realistic lensed images of the 
background elliptical galaxies. The distortions of the lensed images 
have been quantified with the IRAF task ELLIPSE. We find moder- 
ately lensed galaxies can be potentially differentiated from normal 
galaxies as outliers in the coefficients of a^/a^s/a etc. 



24 23 22 21 20 
< fi> e (mag arcsec -2 ) 

Figure 12. Kormendy relation of the elliptical galaxies. The black dashed 
line shows the Kormendy relation of the unlensed background elliptical 
galaxies (black dots). The red dashed line presents the relation of the lensed 
elliptical galaxies (red dots). In this figure, the lens galaxies are at redshift 
Z = 0.3 and the source galaxies at z = 0.6. 



We also explore the offsets in the Kormendy relation and fun- 
damental plane of the lensed elliptical galaxies, and attribute these 
to the magnification effect of the effective radii, while the mean sur- 
face brightnesses within the radii are nearly unaffected. The sys- 
tematic offsets are on the order of 0.03 in logi? e . To observe the 
shift in the fundamental plane, we need to determine the velocity 
dispersion spectroscopically while for the Kormendy relation, we 
only need the source redshift (e.g., from photometric redshift) to 
obtain the physical effective radius, so the observational demand is 
somewhat lower. 

For strong lensing, it is possible to model an individual lens to 
extract information about the lens. For moderate lensing, this will 
be difficult, and its application will be mostly statistical. Since the 
scatter in the fundamental plane is small (0.08 in log 10 W c ), a 0.03 
shift can in principle be deteced with a few thousand galaxies. For 
the Monte Carlo simulations we discussed in ^2] we first fit the un- 
lensed sample with a linear model and then renormalise the error 
bars so that the total x 2 P er degree of freedom is unity. We then fit 
the same model to the lensed sample. If there were no systematic 
shifts, then the resulting x 2 should follow roughly a normal distri- 
bution with mean equal to^ ^ eiin = N& f = N — 2 and a dispersion of 
tr x = y/lNfof, where N is the number of pairs of galaxies. Any sig- 
nificant deviation will indicate the no offset model is unsatisfactory. 
For our simulated sample, we find the x 2 for the no offset model is 
significantly higher than this expectation, and the excess x 2 above 
Xmesm ' s on tne or der of 4.20-f and 6.6cr^ for 1600 and 5000 pairs 
of galaxies. For the Kormendy relation, the corresponding signif- 
icance levels are 'iAcr^ an d 6.3<t m , somewhat lower than that in 
the fundamental plane, which is not surprising given the larger scat- 
ters in this relation. 

In practice, there will be several complications since the lens 
and source redshifts may only be available photometrically. How- 
ever, if there are no systematic errors, we can in principle bin the 
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background galaxies, and stack them to find the systematic offset 
in the fundamental plane and Kormendy relation as a function of 
separation, which in turn provides strong constraints on the aver- 
age profiles of galaxies at large radii. This complements the weak 
galaxy-galaxy method using shear. 

At redshift 0.3, the median Einstein radius is around 0.83 arc- 
sec (corresponding to about 2.6/r 1 kpc), and moderate lensing can 
probe to ~ 50 E , about 13.2/T 1 kpc. This is an interesting radius, 
since it may be close to the regime where the density slope may 
be changing from isothermal (p oc r~ 2 , Koopmans|2009 1 to steeper 
profiles (r~ 3 ): for a galactic-sized halo with virial radius of r v = 200 
kpc and a concentration parameter of c = 10, the radius where the 
density slope changes may be around r v jc = 20 kpc. 

Although the results presented here are promising, there are 
several points which need to be investigated further in the future. 
First, we do not fully consider the effects of seeing, as would 
be necessary in more realistic simulations. However, some of the 
quantities that we use, such as the Fourier components (see Fig.^J, 
are already averaged over a range of radii and so the impact will 
be somewhat limited. Second, while many of the elliptical galax- 
ies can be described well by the de Vaucouleurs profile, others are 
better described by the more general Sersic (1968) profile. Third, 
for simplicity we adopt the singular isothermal lens model. In re- 
ality, several other models can better describe the lens galaxies. 
These include the singular isothermal ellipsoid mode l (|Keeton &| 
|Kochanek|1998) , and the GNFW model ( |Zhao|1996| |Chae|2002> 
which may be particularly suitable for studying the transition in the 
density profiles. In more detailed studies we could parameterise the 
foreground lens using these models and extract the best-fit param- 
eters with more rigourous statistical methods, such as maximum 
likelihood or Bayesian techniques. It may be particular interesting 
to explore what we can learn with photometric redshifts alone for 
moderate gravitational lensing using Monte Carlo simulations. 
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